Geometric chaoticity leads to ordered spectra for randomly interacting fermions 
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A rotationally invariant random interaction ensemble was realized in a single- J fermion model. 
The dominance of ground states with zero and maximum spin was confirmed and explained with a 
statistical approach based on the random coupling of individual angular momenta. The interpreta- 
tion is supported by the structure of the ground state wave functions. 



> 



i 



The interplay of regular and chaotic features in many-body quantum dynamics is currently extensively studied both 
for simple models and for realistic applications to atomic JjJ , nuclear , and condensed matter physics || , as well 
^v^j ■ as for understanding properties of the QCD vacuum Q . Typical finite "shell-model" systems such as complex atoms 
and nuclei are described by the mean field and corresponding residual interaction. The density of the mean field 
configurations grows exponentially for combinatorial reasons, so that the interaction becomes effectively strong at 
sufficiently high excitation energy leading to generic chaotic features both in spectral statistics, which rapidly move to 
the limit of random matrix theory , and in properties of wave functions . Studies of finite many-body systems 
\Q ' have to account for the existence of constants of motion such as total angular momentum, isospin and parity. If these 
conservation laws are exact, one usually deals with the states of each class separately. However, little attention was 
paid to the problem of correlations between classes of states which are described by the same Hamiltonian but belong 
to different values of exact integrals of motion. 

An obvious and practically important example is angular momentum conservation in a finite Fermi-system. The 
7-H , prediagonalization procedure of projecting the correct value J of nuclear spin out of the m-scheme Slater determinants 
induces by itself a strong mixing of the states within a shell model configuration The projected states of various 
spins acquire a nearly uniform degree of complexity and energy dispersion. For a sufficiently large dimension, the 
majority of states correspond to a complicated quasi-random coupling of individual spins. This "geometric chaoticity" 
was used long ago [gj in evaluating the level density for a given J. It also plays an important role in the response 
to external fields, large amplitude collective motion, dissipation and so on The similarity of different J-classes 
with respect to mixing was demonstrated ||,[l0| in the nuclear shell model by the studies of complexity occupation 
numbers, strength functions and pairing properties. This raises also a question of existence of compound rotational 
bands [pr| which would connect complicated states having different J but almost the same mixing. 

A new angle of looking at the problem was introduced by refs. [p^|Jl^ 1 where the spectrum of a random but 
rotationally invariant Hamiltonian was obtained for a shell-model Fermi system. In spite of the random character of 
the two-body interaction, the fraction fo of the ensemble realizations with a ground state spin Jo = was much higher 
than the total statistical fraction /q of J = states in shell-model space. This result was confirmed in refs. |f^ , |f5| , 
as well as for the interacting boson model ]lq] . A new feature discovered in |fj,|l^| was an excess of the probability 
fjmax f° r ground state to have the maximum possible spin J ma x- The emergence of regular features as an output 
of a random interaction seems to contradict the notion of geometrical chaoticity. Below we show that, vice versa, the 
geometric chaoticity provides a base for explaining the main features of the pattern. 

First we give a couple of trivial examples which point out the possible source of the effects, namely an analog of 
the Hund rule in atomic physics. Consider a system of N pairwise interacting spins with the Hamiltonian 

H = AY / Sa-s b = A[S 2 -Ns(s + l)}. (1) 

If the interaction strength A is a random variable with zero mean, then the ground state of the system will have 
equal, fo — fs max = 1/2, probabilities to have spin S = or S = S ma x (antiferromagnetism or ferromagnetism) . A 
similar situation takes place in the degenerate pairing model fl7| where the pair creation, Pq, pair annihilation, Pq, 
and particle number, N, operators form an SU(2) pseudospin algebra. Then the eigenenergy is simply proportional 
to the pairing constant so that, for a random sign of this constant, the ground state pseudospin will be (unpaired 
state of maximum seniority) or maximum possible (fully paired state of zero seniority), on average in 50% of cases. 



In the Elliott SU(3) model 18 1, as well as in any model with a rotational spectrum, the normal or inverted bands will 



happen evenly if the moment of inertia takes positive or negative values randomly. 
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Let us consider a system of interacting fermions. For simplicity we limit ourselves here to a case of N identical 
particles on a single- j shell which provides a generic framework for the extreme limit of strong residual interaction. 
Rotational invariance is preserved, so that all single-particle m-states are degenerate in energy. Within this space, 
the general two-fermion rotationally invariant interaction can be written as 



LA 



where the pair operators with pair spin L and its projection A are defined as 



(2) 



(3) 



and C are the Clebsch-Gordan coefficients. Because of Fermi statistics, only even L values are allowed in the single- j 
space. This fact was ignored in the attempt 12 1 to construct the quasiparticle ensemble with identical distributions of 
the parameters Vl in the particle-particle channel and the parameters Vk for the same interaction transformed to the 
particle-hole channel, H ~ Y]r^ K £k~(q ^a)K K {a 'a)Kii (the difference between the interactions in the two channels was 
discussed long ago by Belyaev |l9| , and served as a justification for an interpolating model "pairing plus multipole- 
multipole forces"). Since K can take both even and odd values, the number of parameters is different in the two 
representations, and Vk cannot be independent if Vl are. 
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b) N=4 j=25/2 
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d) N=6j=17/2 
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f) N=5 bosons 
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FIG. 1. The distribution of ground state angular momenta for various systems of N fermions of spin j, (a-e).The bosonic 
approximation, fj is in panel (J). The dotted lines are the statistical distribution of allowed J and the solid lines are the 
ensemble results. In (c) the dashed line is for Vb = 0, i.e. no pairing. 

Assuming that the coupling constants Vl are random, uncorrelated and uniformly distributed between -1 and 1, we 
get the distribution fj of the ground state spin Jo shown in Fig. [j](a-e) for N = 4 and N — 6 at different values of j. 
For comparison we show by dotted lines the statistical distributions fj based on the fraction of states of given J in 
the entire Hilbert space for given N . The overwhelming probability /o shows the same phenomenon in the uniform 
ensemble as observed earlier in Gaussian ensembles of Vl JlJ,|l^,|lB] . Further evidence of the dominance of Jo = 
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configurations is given by the example, Fig. [I] (e) , for an odd number of particles, where excess of the ground state 
spin Jo = j is evidently related to the ground spin Jo = in the neighboring even system. 

First we note that the effect seems to exist already in a crude approximation modeling fermionic pairs by bosons. 
The commutation relations for the fermion pair operators ([}]) are (L and L' are even), 

[Pl'A' , PL] = 5 L , L 5 A , A + 2 J2 C ^-n C^ala m , . (4) 

mm'n 

The second term in (|]) is of the order N/Cl where Q is the capacity (= 2j + 1 in our case) of the fermionic orbitals. 
It is small for a small number of fermions; for a nearly filled shell its effect is also small because of the particle-hole 
symmetry of states. For intermediate shell occupation this term is not small but can be approximately substituted 
by its mean value (the monopole part with spin K = 0). Then, after a simple renormalization, Pl\ become bosonic 
operators. This is the assumption used in the original boson expansion techniques |20|Jj"l|] and later in the interacting 
boson models: fermionic pairs Pla are substituted by bosons 6la, and the Hamiltonian (0) becomes a sum of random 
bosonic energies X)la The ground state in each realization corresponds to the condensation of the bosons 

into the single-boson states |LA) with the lowest value of oJt- For a given L, the many-boson states with different J 
allowed for the condensate are degenerate, but the value L = is singled out by the obvious fact that for u = min 
all degenerate states have total spin J — while for the minimum boson energy lol at L ^ any specific value of J, 
including J = 0, appears only in a small fraction of states. If all Vl have the same distribution, we expect /o ~ 1/fc 
where A: is a number of (equiprobable) values of L. All other values J ^ appear with small probabilities ~ 1/fc 2 . 
This is demonstrated by Fig. [[](/) where the pattern is qualitatively similar to that in Fig. |j](a-e). The bosonic effect 
gives only a part (decreasing with increasing j) of the Jo = dominance observed for the fermions. Another argument 
against the dominance of the bosonic correlations is given in Fig. 0(c). Here we see that after exact elimination of 
the monopole term (Vl=o = 0), the picture does not significantly change although the value Vo is now the lowest only 
in a small fraction, ~ 2~( fc_1 ), of all cases (when all Vl^o are positive). 

In our opinion, the main effect comes from the statistical correlations of the fermions. They resolve the bosonic 
degeneracy in favor of the J = and J = J max ground states. In the strong mixing among nearly degenerate states, 
the eigenstates emerge as complicated chaotic superpositions. The only constraints left are the conservation laws for 
the particle number and total spin. The latter can be taken into account by the standard cranking approach |j|, p2"| , p3"t . 
Thus, we model the system by the Fermi-gas in statistical equilibrium with the occupation numbers n m of individual 
orbitals characterized by the angular momentum projection m onto the cranking axis. The presence of the constraints 
creates a "body-fixed frame" and splits effective quasiparticle energies, although instead of the collective rotation 
around a perpendicular axis we have here a random coupling of individual spins with the symmetry (cranking) axis 
being the only direction which is singled out in the system |24[ . Under the constraints 

N = ^n m , M = y^mn mi (5) 



equilibrium statistical mechanics leads to the Fermi-Dirac distribution 

1 



exp(7m — (j,) + 1 



(6) 



determined by the Lagrange multipliers of the chemical potential /z and cranking frequency 7; in the end the total 
projection M (equivalent to the K quantum number for axially deformed nuclei) is identified with the total spin J. 

The quantities fi(N, M) and 7 (AT, M) can be found directly from (||). At M = we have 7 = 0, so that the expansion 
in powers of 7 allows one to study the most important region around M — 0; the power expansion is sufficient for all 
M except for the edges. With no cranking, one has the uniform distribution of occupancies = n = N/Cl. With 
the perturbational cranking, the occupation numbers are 



7 2 

1 - 7771(1 - n) + y (™ 2 - (fn 2 )){l - n)(l - 2n) + 



(J) 



Here (m 2 ) = (l/^)X) m m2 = J 2 /3, and terms of higher orders are not shown explicitly. The expectation value of 
energy in our statistical system can be written as 

(H)= Yl V L \C% m f(n mi n m2 ). (8) 

LArnim2 
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Neglecting the correlations between the occupation numbers, (n rni n m2 ) w n mi n m2 , we use the statistical result (ffl) 
and calculate the geometrical sums with the Clebsch-Gordan coefficients. Expressing the parameter 7 in terms of the 
total spin M — > J, we come to the result including the terms of the second order in J 2 , 

{H) N .j = ^(2L + l)V L [h (L) + h 2 (L)J 2 + hi{L)j\ (9) 

L 

where 

h (L)=n 2 , ft a (L) = |^^ ) (10) 



9 (1 - 2n) 2 (3L 4 + 3L 2 - 12j 2 L 2 - 6j 2 + 8j 4 ) 
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(11) 



Jo is determined by the ensemble distributions of h 2 ,4 = ^2l(2L + tyVchz 4,(L) . For all realizations of the random 
interaction with non-negative h 2 and positive /14, the ground state has spin J Q — 0. If /12 > but /14 < 0, one 
has a local minimum of energy at J — although there is a possibility to reach the absolute energy minimum at 
Jmax = (l/2)JV(f2 — -/V). This will not happen if at J = J maa; we still have h 2 + Jmax^A > 0- Therefore the probability 
to have the ground spin state equal to zero turns out to be, in this approximation, 



S(h 2 .hi) 



dh 2 dhiP '{h2hi) , 



(12) 



where the region S is defined by the conditions h 2 > 0, /i 4 > —(h 2 / 'Jmax)- Since /i 4 is small, the result is close to that 
for the right semi-plane h 2 > 0, and fo should be close to 50%. For a Gaussian distribution of the parameters Vl 
with zero mean and variances <7l, the distribution of the linear combinations /i2,4 is again Gaussian, and the integral 
over the region S in (|l2|) gives for this case 



fo = 7 + 7^-arctan 

4 Z7T 
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(13) 



which is close to 1/2. Here we introduced the combinations of geometric factors weighted with the corresponding 
variances, 



A = 5> 2 (L))V 2 , 



D = J2h 2 (L)h i (L)a 2 , 



B = Y J {h i {L)) 2 a 2 L . 



(14) 



The 7-expansion fails for large momenta. However, the states with high M can be constructed exactly. For Fig. 
H we used our statistical approach near J = in conjunction with the exact values in the end region J = J max 
to improve the above result for fo and to get an upper bound for /j max . Thus the statistical approach provides a 
good estimate for the dominance of J = and J = Jmax in the ground state; more subtle effects such as odd-even 
staggering should be considered separately. 
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FIG. 2. fo for N — 4 and different j; ensemble results (solid line), statistical theory (dotted line, panel (a); upper limit for 
fj m ai from the statistical theory and analysis of the Jmax region (dotted line, panel (b)). 
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FIG. 3. The distribution of overlaps of Jo = ground states of the degenerate pairing model (Vb = —1, Vl^o = 0) with 
those for the ensemble choices a) random ensemble with Vb = — 1, b) random ensemble. N = 6 and j — 11/2 in both cases, 
and the dotted line is the predicted P(x). 

Although the energy spectra with random two-body interactions bear clear resemble the ordered spectra of pairing 
forces, the structure of the eigenstates is close to that expected for chaotic dynamics jL4|. Fig. ||(b) shows the 
distribution P(x) of the overlaps x = |(J — 0, g.s.|0, p) | 2 of ground states with spin obtained in the random 
ensemble with the ground state |0,p) for the degenerate pairing model, the latter corresponding to the case of fixed 
Vq = — X, Vl^o = 0. In the chaotic limit the wave functions are expected [0,^2j to behave as random superpositions 
of basis states with uncorrelated components C uniformly spread over a unit sphere, P(C) cx 8(^2 C 2 — 1). This is 
equivalent to the distribution of a single component P{C\) <x (1 — C 2 )( n ~ 3 ^ 2 where n is the space dimension. For 
n>l, the distribution P{C\) is close to Gaussian whereas the overlaps x — Cf obey the Porter-Thomas distribution. 
In the case of Fig. || (N = 6 particles, j = 11/2) the dimension of the J = space is small, n — 3, so that P{C\) 
is constant, and we expect P(x) cx as in the case of the pion multiplicity for the disordered chiral condensate. 

Another case considered in Fig. ^|(a) corresponds to the overlap of the degenerate pairing model ground state with 
the ground state in the model with Vq = —1,Vl=£o random. Of course, here the completely paired state can appear as 
the ground state even for random strengths in the channels L =/= which gives the peak at the overlap x = 1. But the 
character of the distribution changes as well becoming effectively two-dimensional: for n — 2, P{x) cx 1/ *Jx(l — x). 

To conclude, we have shown that statistical correlations of fermions in a finite Fermi system with random inter- 
actions drive the ground state spin to its minimum or maximum value. This effect is related to the geometrical 
chaoticity of the random spin coupling of individual particles. This means, that the dominance of + ground states 
in even-even nuclei may at least partly come from incoherent interactions rather than solely from coherent pairing. 
The structure of ground states with an "antiferromagnetic" type ordering, Jo = 0, is compatible with the predictions 
for chaotic dynamics. Quantitative relations between the effects of geometric chaoticity and pure dynamic effects in 
finite many-body systems should be an interesting subject for further detailed studies. 
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